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j  13.  ABSTRACT  (Mdjamum  200  wort*)  1 

It  has  been  demonstrated  that  in  a  sooting  laboratory  flame  the  time  to  add  a  carbon  atom  to  a  growing  carbon 
species  is  about  the  same,  or  shorter,  for  the  ionic  than  for  the  free  radical  mechanism.  The  calculated  step 
time  is  consistent  with  the  experimentally  observed  time  in  the  flame  to  build  a  soot  particle.  It  has  also  been 
demonstrated  that,  up  to  about  30  carbon  atoms,  the  thermodynamic  driving  force  for  the  ionic  mechanism 
is  much  greater  than  for  the  neutral  mechanism;  above  30  carbon  atoms  they  are  equal.  The  Langevin  theory 
of  ion-molecule  reactions  has  been  modified  to  accommodate  large  ions.  The  results  indicate  a  greater 
collision  coefficient  for  large  ions  than  calculated  by  Langevin.  The  coefficient  for  a  hard  sphere  collision 
decreases  as  ion  diameter  increases  up  to  about  3  nm  diam  for  a  conducting  spherical  ion,  and  about  1  nra 
for  a  nonconducting  spherical  ion.  The  Sandia  Flame  Code  has  been  modified  to:  (1)  accept  individual 
experimental  concentration  profiles  as  input;  (2)  use  non-Arrhenius  rate  constants;  and  (3)  accept  ambipolar 
diffusion  coefficients.  Analysis  of  our  previous  modeling  of  the  ionic  mechanism  and  many  runs  with  the 
modified  code  indicate  either  a  basic  flaw  in  the  mechanism  or  a  modeling  problem.  The  detailed  modeling 
results  are  inconsistent  with  the  calculated  speed  with  which  carbon  atoms  are  added  and  with  the  strong 
thermodynamic  driving  force  for  the  ionic  mechanism. 
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I.  INTRODUCTION 


This  is  the  first  annual  report  on  a  program  to  compare  the  ionic  mechanism  of  soot  formation 
with  a  free  radical  mechanism  using  computer  modeling.  In  a  previous  contract1  with  the  same  objective 
we  collaborated  with  Michael  Frenklach  and  Hai  Wang  at  Penn  State.2  In  that  program  we  developed 
the  detailed  chemical  kinetic  scheme,  estimated  the  elementary  reaction  rates  involved,  and  developed  the 
thermodynamics  of  the  ions  involved  in  the  ionic  mechanism;  this  information  had  been  unavailable. 
Frenklach  and  Wang  improved  the  free  radical  mechanism  and  incorporated  our  ionic  mechanism  into 
the  Sandia  Flame  Code  which  they  used  for  the  free  radical  mechanism.  They  ran  the  program  on  an 
IBM  3090  computer. 


They  experienced  several  problems  in  adapting  the  Sandia  code  to  handle  ions:  (1)  ionic  diffusion, 
usually  treated  as  ambipolar  diffusion  could  not  be  incorporated  into  the  code  so  an  artificial  fit  was 
employed-this  was  accurate  enough  for  these  calculations,  but  was  conceptually  unpleasant;  (2)  non- 
Arrhenius  ion-molecule  reaction  rate  coefficients  could  not  be  accommodated  by  the  code;  and  (3) 
experimental  ion  and  neutral  profiles  could  not  be  used  as  input  to  the  code.  The  third  limitation  was 
the  most  serious  because  the  program  calculated  excessive  concentrations  of  small  neutral  species  involved 
in  the  chemiionization  steps  and  thus  calculated  too  high  concentrations  of  initial  ions.  It  also  had  to 
carry  along  a  large  number  of  neutral  species  which  made  no  contribution  to  the  ionic  mechanism  and 
thus  wasted  the  computer  resources.  These  things  limited  our  capability  of  evaluating  the  ionic 
mechanism,  and  the  code  predicted  greater  ion  concentrations  than  would  have  been  obtained  were  the 
neutral  concentrations  consistent  with  experimental  data.  Another  major  problem  involved  turnaround 
time;  it  usually  took  weeks  to  get  a  revised  mechanism  run,  thus  greatly  limiting  the  number  of  computer 
experiments  that  could  be  performed. 


To  improve  the  situation,  we  hired  Fokian  Egolfopoulos  (originally  at  Princeton  University,  now 
at  the  University  of  Southern  California)  to  adapt  the  Sandia  code  to  ions.  He  was  very  successful  in  this 
endeavor;  the  three  above-mentioned  problems  were  corrected.  AeroChem  purchased  a  Silicon  Graphics 
workstation  on  which  it  installed  the  revised  flame  code.  This  operates  very  satisfactorily,  with  run  times 
typically  about  15  min  to  two  hours  and  turnaround  time  instantaneous.  This  has  allowed  us  to  make 
more  computer  runs,  whereby  we  have  identified  problems  with  the  mechanism,  or  with  the  code,  which 
were  previously  unsuspected. 

During  this  report  period  we  have,  in  addition  to  the  above:  (1)  continued  to  expand  the 
thermodynamic  data  base  on  ions3;  (2)  completed  modifications  to  the  Langevin  theory  of  ion-molecule 
reactions  to  accommodate  large  ions  (contrary  to  expectations  the  rate  coefficients  were  increased);  (3) 
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fine  tuned  the  comparison  between  the  free  radical  mechanism  and  the  ionic  mechanism  by  calculating 
the  times  to  add  a  carbon  atom  to  the  growing  species,' 4,5  (the  ionic  mechanism  appears  to  be  faster);  (4) 
compared  the  free  radical  mechanism  with  the  ionic  mechanism  by  determining  the  thermodynamic 
driving  force  between  the  two  mechanisms  (the  driving  force  for  the  ionic  mechanism  is  much  greater); 
and  (5)  organized  the  experimental  measurements  made  on  benzene/oxygen  flames  under  a  previous 
contract6  to  form  the  experimental  data  base  for  comparing  a  computer  model  of  the  benzene/oxygen 
flame  (this  flame  is  much  different  from  the  acetylene/oxygen  flame). 


fATEMENT 


The  ionic  mechanism  of  soot  formation  in  flames  will  be  further  evaluated  and  compared  with 
the  neutral  mechanism  by  pursuing  the  following  phases: 

Phase  I.  Extend  the  Ionic  Mechanism  to  Benzene-Oxygen  Flames 

1 .  Extend  the  thermodynamic  and  reaction  rate  coefficient  data  base  to  include  those 
species  found  in  benzene/oxygen  flames  which  are  not  present  in 
acetylene/oxygen  flames. 

2.  Organize  the  experimental  data  available  on  the  benzene  flame  including  both 
neutral  and  ionic  species. 

3.  Develop  a  detailed  ionic  mechanism  for  the  formation  of  soot  in  benzene  flames 
and  submit  this  to  others  to  run  on  a  large  computer  to  compare  its  agreement 
with  experimental  data  and  to  compare  its  simulation  of  soot  formation  with  the 
neutral  mechanism. 

4.  Compare  the  computer  modeling  results  with  experimental  data  and  interpret  the 
results  in  terms  of  the  major  chemical  pathway  to  soot  in  the  benzene  flame. 

Phase  II.  Model  Coagulation  and  Agglomeration 

1 .  Extend  the  detailed  mechanism  of  soot  nucleation  to  include  coagulation  of  large 
ions  and  neutrals,  and  charged  with  charged  and  neutral  incipient  particles. 
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2.  Program  a  desktop  computer  to  test  the  model  developed  in  Task  1  (in  a  limited 
fashion)  against  experimental  data.  If  warranted,  submit  this  model  to  others  to 
test  on  a  mainframe  computer. 

3.  Interpret  the  results  from  Task  2  in  terms  of  a  simplified  mechanism. 

Phase  HI.  Develop  a  Theory  for  Large  Ion-Molecule  Reactions 

1.  Extend  the  Langevin  theory  of  ion-molecule  reactions  to  include  large  ions  by 
removing  the  restriction  of  a  point  charge  on  the  ion. 

Phase  IV.  Compare  Thermodynamic  Predictions  of  Soot  Formation  with  Experimental 
Observations 


1 .  Collect  and  organize  the  literature  data  on  soot  yield  and  acetylene  concentrations 
as  a  function  of  equivalence  ratio. 

2.  Calculate  the  thermodynamic  equilibrium  concentration  of  acetylene  and  soot  as 
a  function  of  equivalence  ratio  for  the  literature  systems  identified  in  Task  1. 

3.  Compare  the  experimental  and  calculated  values  and  interpret  them  in  terms  of 
generalizations  relevant  to  the  mechanisms  of  soot  formation. 


HI.  THERMODYNAMICS 


Thermodynamic  data  continue  to  be  developed  for  ions  involved  in  the  acetylene/oxygen  flame. 
A  paper3  containing  die  thermodynamic  data  we  have  developed  to  date  has  been  submitted  to  the  Journal 
of  Physical  and  Chemical  Reference  Data  for  publication.  In  the  near  future  we  will  initiate  collection 
and  development  of  data  on  the  oxygenated  ions  that  appear  in  benzene  flames  but  not  in  acetylene 
flames.  This  will  be  done  in  preparation  for  modeling  the  benzene/oxygen  flame. 
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IV.  REACTION  MECHANISM  AND  REACTION  COEFFICIENTS 


The  reaction  mechanism  and  reaction  rate  coefficients  previously  reported1  have  been  changed 
very  little. 


V.  iAK.GEVIN.THlQllY.  ..EQBJLABGE  IQN§ 


We  have  used  the  Langevin  theory7  to  calculate  ion-molecule  reaction  rate  coefficients.  This 
theory  has  repeatedly  been  shown  to  agree  with  experiments  for  small  ions  near  room  temperature,  but 
has  not  been  tested  for  large  ions  at  flame  temperatures.  The  theory  assumes  a  point  electric  charge, 
which  for  large  ions  seems  questionable.  In  Langevin’s  theory  the  collision  rate  constant  is  independent 
of  molecular  velocity,  and  thus  temperature.  However,  more  conventional  neutral  collision  rates  increase 
with  temperature.  The  objective  of  the  present  study  is  to  generalize  the  Langevin  theory  to  account  for 
ions  of  finite  dimensions.  The  extension  to  finite  dimensions  results  in  a  temperature  dependence.  A 
temperature  effect  has  been  investigated  by  Bei  et  al.8  and  others  for  small  ions  reacting  with  permanent 
dipole  molecules,  but  the  effect  of  temperature  has  not  been  demonstrated  for  nonpolar  molecules  and 
has  not  been  studied  for  large  ions. 

Our  work  is  being  prepared  for  publication  so  it  will  be  summarized  here  and  some  results 

given.* 


Three  ion  models  were  studied.  The  first  model  considers  the  electronic  charge  to  be  at  a  point 
at  the  center  of  a  nonconducting  sphere  of  finite  radius  whose  dielectric  constant  equals  the  vacuum  value. 
Thus  the  electric  field  distribution  is  identical  to  Langevin’s,  but  additional  collisions  are  possible  for 
faster  relative  speeds  due  to  the  physical  size  of  the  ion.  This  model  is  similar  to  that  of  Ref.  8. 

The  second  model  considers  the  ion  to  be  a  conducting  sphere  of  finite  radius.  For  this  ion  in 
free  space  the  charge  is  uniformly  distributed  over  the  surface  of  the  sphere.  When  polarizable  neutral 
molecules  are  present,  the  neutrals  become  dipoles  under  the  influence  of  the  ion  electric  field,  and  they 


*See  Section  IV  for  Nomenclature. 
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act  to  redistribute  the  ion’s  surface  charge.  This  redistribution  of  the  ion  surface  charge  increases  both 
the  electric  field  strength  at  the  neutral  and  its  induced  dipole  moment. 

The  third  model,  "ion  polarizability,"  considers  the  ion  to  be  a  nonconducting  sphere  of  finite 
radius  with  the  ionic  charge  in  the  center  as  in  the  first  case,  but  with  the  capability  of  supporting  an 
induced  dipole  moment  at  the  sphere’s  center  under  the  influence  of  an  external  electric  field.  A  fourth 
approach  would  be  to  displace  the  charges  making  up  the  induced  dipole  within  the  sphere,  but 
insufficient  information  is  available  to  determine  the  displacement. 


A.  NONCONDUCTING  SPHERE 

A  trajectory  analysis  using  standard  relationships  for  the  electrostatic  force  on  an  induced 
molecular  dipole  leads  to  the  following  relation  between  the  impact  parameter  b  and  the  minimum 
separation  distance  between  molecules  that  do  not  collide 


b 

*Vnin 


2 


j  +  ?£(rmin) 


(1) 


In  Langevin’s  theory  a  critical  radius  rc  exists  for  which  r  >  rc  corresponds  to  trajectories  that 
do  not  lead  to  particle  capture  and  r  <  rc  corresponds  to  capture  with  r  always  decreasing  monotonically 
to  zero.  The  limiting  case  r  =  rc  corresponds  to  orbital  motion.  In  this  case  the  electrostatic  and 
centrifugal  forces  are  balanced  so  that 


firi)2 


(2) 


Using  the  conservation  of  energy  and  angular  momentum  for  the  molecular  trajectory  and 
employing  Eq.  (2)  leads  to  the  following  expression  for  rc: 
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r  pdB  (3) 

^  dr  2  y 

Equation  (3)  is  the  cornerstone  of  our  analysis  independent  of  the  ionic  model  chosen.  The 
assumption  is  that  the  induced  dipole  axis  is  always  aligned  with  r,  the  vector  connecting  the  centers  of 
the  ion  and  the  neutral  molecule,  throughout  the  trajectory  from  infinite  radius  to  rc.  In  our  numerical 
computations  the  mean  value  of  the  three  principal  polarizability  components  for  acetylene  was  used.  The 
ratio  of  the  mean  to  the  minimum  value  is  1.187.9 

Once  an  ion  model  is  chosen  one  can  calculate  E  and  dE/dr  as  functions  of  r,  and  re  is  that  value 
of  r  that  satisfies  Eq.  (3)  for  a  given  value  of  va  and  neutral  molecular  properties.  Setting  rc  =  r^  in 
Eq.  (1)  results  in  a  value  for  the  impact  parameter  b  in  terms  of  the  relative  system  kinetic  energy  at 
infinite  separation.  All  trajectories  with  b  less  than  this  value  will  lead  to  collisions.  Finally,  the  total 
collision  rate  k  can  be  determined  from 


k 


(  L  2  V 


(4) 


with  the  Maxwell-Boltzmann  distribution  in  kinetic  energy  used  for  dn^/n. 

This  treatment  is  applicable  to  a  nonconducting  sphere.  To  this  point  there  is  no  real  difference 
between  this  formulation  and  Langevin’s  if  Eq.  (3)  alone  is  used  to  find  the  impact  parameter.  However, 
Eq.  (3)  is  not  relevant  when  rc  <  a;  +  a,,,  where  a;  and  a„  are,  respectively,  the  ion  and  neutral  molecule 
radii.  In  this  case  collisions  occur  at  a  given  value  of  v*  for  all  impact  parameters  b  less  than  that  given 
by  Eq.  (1)  with  r^  =  %  +  a„. 

Equation  (4)  is  integrated  in  two  kinetic  energy  ranges  separated  by  the  condition  rc  =  aj  +  a„. 
For  the  upper  range  of  rc  the  Langevin  mechanism  holds  independent  of  the  ion  or  neutral  radii,  but  for 
smaller  values  of  re  there  are  additional  collisions  not  included  in  the  Langevin  theory.  This  additional 
class  of  collisions  becomes  more  important  with  increasing  temperature.  Even  in  this  case  the  trajectories 
are  still  affected  by  the  ion  electric  field  so  that  the  impact  parameter  is  larger  than  for  strictly  neutral- 
neutral  collisions. 
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B.  CONDUCTING  SPHERE 

For  the  second  ion  model  of  a  conducting  sphere  we  determine  the  effect  of  a  dipole  on  the 
sphere’s  surface  charge  distribution  in  terms  of  the  simpler  solution  for  an  external  point  charge  on  a 
charged  conducting  sphere. 

The  classic  problem10  is  that  of  a  point  charge  of  strength  qp  external  to  a  conducting  sphere  of 
radius  %  and  total  integrated  surface  charge  the  ionic  charge.  In  the  absence  of  an  external  charge 
the  electric  field  can  be  represented  by  a  point  charge  of  strength  %  at  the  center  of  the  sphere.  With  an 
external  charge  present  an  image  method  is  used  to  construct  the  electric  voltage  potential  field  exterior 
to  the  sphere  in  terms  of  a  simple  distribution  of  point  sources  in  the  sphere.  The  external  source  induces 
an  image  source  within  the  sphere  located  along  a  line  joining  the  center  of  the  sphere  and  the  external 
point  charge  at  a  distance  af/y  from  the  center,  where  y  is  the  distance  from  the  center  of  the  sphere  to 
the  external  point  charge.  The  charge  of  the  image  source  is  \  qp/y.  To  maintain  the  net  ionic  charge 
the  value  of  the  image  charge  is  subtracted  from  the  value  of  at  the  sphere’s  center. 

The  electric  field  due  to  a  dipole  and  a  conducting  sphere  is  obtained  by  superimposing  the  above 
solution  for  a  positive  point  charge  with  that  of  a  negative  charge,  displaced  from  the  positive  charge  to 
form  the  dipole.  The  resulting  field  is  quite  different  from  that  in  Langevin’s  model  when  the  dipole  and 
sphere  are  close  to  each  other.  The  field  approaches  infinity  at  a  separation  distance  r,  which  can  be 
greater  than  the  sum  of  the  ion  and  molecule  radii.  With  an  infinite  field  a  collision  occurs  independent 
of  kinetic  energy  so  that  the  equivalent  hard  sphere  radius  of  the  conducting  spherical  ion  is  greater  than 
its  physical  radius. 

In  our  computations  Eqs.  (3)  and  (1)  and  the  new  expressions  for  E  are  used  with  r^  =  rc  when 
rc  >  r,.  When  rc  <  r„  Eq.  (1)  is  used  to  determine  the  impact  parameter  with  r^  =  rf. 

C.  ION  PQLAMZASIUTY 

The  third  model  considers  the  ion  as  well  as  the  approaching  neutral  to  be  polarizable.  Thus, 
the  charge  on  the  ion  induces  a  dipole  moment  on  the  approaching  neutral  which  in  turn  induces  an 
electric  field  on  the  ion,  inducing  a  dipole  moment  on  the  ion. 

The  polarizability  of  aromatic  hydrocarbon  ions  has  been  studied11  using  molecular  orbital  theory 
and  compared  to  data  for  naphthalene  (C10Hg).  In  this  case  polarizability  for  positive  and  negative  ions 
is,  respectively,  within  about  4  and  18%  of  the  value  for  the  neutral  species.  For  our  present  study  the 
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ton  dipole  will  be  assumed  to  be  located  at  the  center  of  the  ion.  It  would  be  more  realistic  to  have  a 
charge  distribution  on  the  ion  surface  consistent  with  the  dipole  strength,  but  we  have  not  attempted  to 
model  this. 


Appropriate  treatment  leads  to: 


*20  - 


(5) 


and  this  and  its  derivative  are  used  in  Eq.  (3)  to  determine  the  impact  parameter.  Although  it  is  possible 
for  the  denominator  of  Eq.  (5)  to  be  singular,  we  have  not  found  this  to  be  the  case  for  the  range  of 
molecular  parameters  of  interest. 

With  the  electric  field  specified,  the  impact  parameter  is  determined  by  Eqs.  (1)  and  (3)  when 
rc  >  a*  +  a„  with  rmin  =re,  and  by  Eq.  (1)  for  rc  <  a*  +  a„  with  r^  =  a;  +  a„.  Integration  over  the 
Maxwell-Boltzmann  distribution  is  split  into  two  regions  of  velocity  space  to  reflect  the  two  different 
models  of  impact  parameter. 


D.  RESULTS  AND  PISCUSSIQN 

Results  are  shown  in  Fig.  1  for  the  ratio  of  the  expanded  Langevin  collision  rate,  kEL,  calculated 
as  above,  to  the  Langevin  collision  rate,  kL,  as  a  function  of  temperature  for  the  conducting  sphere  and 
nonconducting  sphere  cases.  The  neutral  molecule  is  assumed  to  be  acetylene,  and  the  ion  is  considered 
to  be  a  sphere  with  a  density  of  1.5  g/cm3.  Curves  are  shown  for  ion  diameters  of  0.5,  1,  and  2  nm. 

It  is  seen  in  Fig.  I  that  the  collision  rate  for  the  expanded  Langevin  theory  is  increased  by  about 
five  over  die  Langevin  rate  for  the  2  nm  nonconducting  ion  at  typical  flame  temperatures.  The 
conducting  sphere  case  leads  to  an  additional  increase  by  nearly  a  factor  of  2.  The  general  increase  with 
temperature  is  due  to  the  finite  size  of  the  ions  presenting  a  hard  sphere  cross-section. 

This  latter  effect  is  illustrated  in  Fig.  2  in  which  the  ratio  of  calculated  collision  rate  to  hard 
sphere  collision  rate  is  plotted  as  a  function  of  temperature.  There  is  little  dependence  on  temperature 
at  higher  temperatures;  the  temperature  effect  is  less  for  the  larger  ions.  In  particular,  the  2  nm 
nonconducting  case  is  very  close  to  the  hard  sphere  neutral-neutral  collision  value  at  all  temperatures. 
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The  effect  of  the  ion  diameter  on  the  ratio  of  the  expanded  Langevin  rate  to  the  hard  sphere  rate 
is  shown  in  Fig.  3.  As  the  diameter  increases,  the  enhanced  rate  coefficient  due  to  an  electric  charge 
becomes  less. 

The  third  case,  of  an  induced  dipole  in  the  ion,  is  not  shown  because  the  data  differed  little  from 
the  nonconducting  sphere  case,  even  for  the  lowest  temperatures  considered.  The  reason  is  that  the  dipole 
is  assumed  to  be  located  at  the  sphere’s  center  so  that  the  separated  ionic  dipole  charges  cannot  get  very 
close  to  the  neutral  molecule  dipole.  Note  that  the  x~6  dependence  in  Eq.  (5)  typical  of  dipole-dipole 
interactions,  means  that  the  effect  falls  off  very  rapidly  with  distance.  If  this  case  were  modeled  with 
a  finite  charge  separation,  the  effect  would  be  much  larger. 

In  this  study  we  have  looked  at  generalizations  of  the  Langevin  theory.  An  important 
consideration  not  treated  is  that  of  geometry  since  large  polycyclic  aromatic  hydrocarbons  are  often  in 
the  form  of  disks  rather  than  spheres  as  assumed  here.  For  example,  the  collision  rate  for  a  disk  having 
the  same  mass  as  a  sphere  can  be  much  greater  than  the  sphere’s  rate  due  to  the  increase  in  surface  area. 
This  is  easily  shown  by  taking  the  rate  k,  for  a  sphere  to  be: 


*a,vm 


(6) 


where  a,  is  the  sphere  radius  and  vm  is  the  mean  molecular  speed.  For  a  disk  of  radius  ad,  considering 
only  collisions  with  the  flat  sides  and  not  the  edge,  the  sum  of  the  collision  rates  kd  for  both  sides  of  the 
disk  is: 


(7) 


If  the  disk  has  a  thickness  h,  then  equating  the  volumes  of  the  sphere  and  disk  leads  to 


k_d_  _  H 
'  3  h 


or 
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h 

(9) 


Extending  the  analysis  to  compute  the  electric  field  of  nonspherical  ions,  e.g.,  conducting  and 
nonconducting  prolate  spheroids  and  disks  is  possible. 12  However,  the  trajectory  analysis  would  be  much 
more  complex. 

Other  considerations  are  to  use  a  fixed  dielectric  constant  for  the  ion,  which  should  yield  results 
intermediate  to  those  of  the  conducting  and  nonconducting  cases,  and  to  consider  an  off  center  location 
for  the  net  ionic  charge.  These  ideas  can  be  extended  to  the  nonspherical  cases  as  well. 

In  summary,  we  have  shown  that  treatment  of  finite  size  ions  leads  to  collision  rates  in  excess 
of  those  predicted  by  classical  Langevin  theory.  The  direction  of  future  work  in  this  area  will  be  guided 
by  our  ability  to  model  the  significant  structural  and  electrical  properties  of  large  hydrocarbon  ions.  This 
analysis  has  not  yet  been  incorporated  into  the  mechanism. 


VI.  COMPUTER  MODEL 


In  our  previous  work  our  model  was  run  on  an  IBM  3090  computer  by  M.  Frenklach  and  Hai 
Wang  at  Penn  State.2  There  was  a  long  turnaround  time  and  their  code  was  not  adaptable  to  ions  or 
some  of  the  computer  tests  which  we  wanted  to  do.  We  thus  hired  Fokian  Egolfopoulos  at  Princeton 
University  (now  at  University  of  Southern  California)  as  a  consultant  to  modify  the  flame  code  package, 
based  on  Sandia’s  CHEMKIN  and  Flame  Codes.  Modifications  were  made  in  the  code  to  incorporate: 
(1)  ion  chemistry;  (2)  ambipolar  diffusion  of  ions;  (3)  the  ability  to  input  neutral  and  ion  profiles;  and 
(4)  the  ability  to  employ  a  modified  Arrhenius  form  for  entering  forward  rate  coefficients. 

The  modifications  to  include  ambipolar  diffusion  coefficients  replace  the  previous  technique  of 
using  neutral  diffusion  coefficients  of  similar  sized  molecules.  Brown  in  his  treatment,1  used  Poisson’s 
equation,  which  is  the  more  correct  way  to  account  for  diffusion.  In  addition,  the  revised  code 
incorporates  a  charge  neutrality  check.  Providing  for  incorporation  of  ion  or  neutral  species  profiles  in 
the  code  allows  us  to  employ  as  input  experimental  ion  concentration  profiles  for  several  initial  species 
(we  have  thus  far  used  only  one  ion),  and  experimental  concentration  profiles  for  neutral  species  profiles 
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important  in  the  ionic  mechanism.  This  is  similar  to  the  frequently  used  technique  of  inputing 
experimental  temperature  profiles  as  data.  It  overcomes  two  problems  previously  experienced  when 
evaluating  the  ionic  mechanism.  First  we  had  to  include  the  neutral  species  mechanism  for  small  species 
to  account  for  the  production  of  chemiions.  Thus  we  originally  included  about  40  neutral  species  and 
about  200  reactions  from  the  neutral  mechanism  of  soot  formation.  This  added  to  the  complexity  of  the 
program  and  increased  computation  time.  The  second  problem  was  that  the  neutral  mechanism  produced 
an  excess  of  chemiions  over  the  experimentally  observed  quantities.  We  were  thus  biasing  the  ionic 
mechanism  and  were  not  testing  the  ionic  mechanism  alone  but  rather  a  combination  of  the  ionic 
mechanism  and  the  neutral  mechanism.  By  using  experimental  data  as  input  we  can  more  accurately  test 
the  ionic  mechanism. 

In  the  present  model,  the  ion-molecule  reaction  rate  coefficients  in  the  forward  direction  are  fixed 
so  that  the  equilibrium  constant  of  the  reverse  reaction  never  exceeds  the  Langevin  rate.  When  the  free 
energy  of  the  reaction  changes  sign  through  the  temperature  range  of  importance  this  is  not 
straightforward.  Therefore  a  modified  Arrhenius  rate  coefficient  equation  was  added  to  the  flame  code 
to  accomplish  this.  This  was  done  by  making  the  forward  rate  coefficients  for  those  reactions: 


k  -  A  x  Tn  x  exp  (—EactJRT)  +  B 


(10) 


The  constants:  A,  B,  Ewt,  and  B  are  determined  empirically  by  fitting  Eq.  (10)  to  a  curve  of  the  forward 
rate  plotted  against  the  temperature,  maintaining  either  the  forward  or  reverse  rate  at  the  Langevin  rate 
via  the  equilibrium  constant.  The  computer  code  had  to  be  modified  to  accept  Eq.  (10). 

To  test  the  computer  models,  AeroChem  purchased  a  Silicon  Graphics  IRIS  workstation.  Run 
times  have  been  relatively  short,  from  a  few  minutes  to  several  hours. 


VII.  COMPARISON  OF  CALCULATED  AND  EXPERIMENTAL  ION 
CONCENTRATION  PROFILES 


The  availability  of  the  modified  Sandia  flame  code  to  handle  ions  and  the  availability  of  a 
workstation  have  greatly  accelerated  the  rate  of  progress  on  this  program.  As  already  implied,  some  of 
this  seems  to  be  in  the  wrong  direction.  There  are  several  problems.  Thermodynamic  equilibrium  seems 
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to  have  a  much  greater  effect  on  the  ion-molecule  reactions  than  expected  (the  reverse  reactions  are 
calculated  from  the  forward  rate  and  the  equilibrium  constant);  this  shows  up  as  a  stronger  temperature 
dependence  than  anticipated.  At  higher  temperatures  the  ion-molecule  reverse  reactions  become  more 
important  because  entropy  effects  become  more  important  with  increasing  temperature: 


A  Gf  =  AH f-  TASf 


(11) 


This  same  effect,  of  course,  controls  neutral  mechanisms  and  may  very  well  explain  why  the  rate 
of  soot  formation  peaks  in  the  temperature  range  of  1400  to  1900  K  (see,  e.g..  Ref.  13).  The  usual 
explanation  involves  a  competition  between  the  rate  of  soot  formation  and  the  rate  of  soot  oxidation. 

The  abnormal  temperature  effect  on  ion  formation  has  shown  up  in  several  ways.  In  some  of  our 
first  computer  runs  with  Frenklach  and  Wang,  a  table  with  the  experimental  flame  temperature  vs. 
distance,  in  which  one  point  was  in  error,  was  inadvertently  used  as  input.  This  showed  up  in  the 
calculation  of  ion  profiles  as  a  dip  in  the  concentration  profiles.  In  the  ion  profile  calculations  made  by 
Frenklach  and  Wang,  the  ions  always  peaked  much  closer  to  the  burner  than  the  experimental  data,  see 
e.g..  Fig  4;  however,  the  maximum  calculated  concentrations  were  very  close  to  the  experimental 
maxima.  We  now  interpret  this  as  a  combination  of  the  temperature  effect  (the  temperature  is  lower  in 
the  early  part  of  the  flame)  and  the  use  of  calculated  rather  than  experimental  neutral  species  profiles  in 
the  ionic  mechanism.  We  currently  use  experimental  concentration  profiles  for  those  species  which  enter 
into  the  ionic  mechanism.  Using  Brown’s  code,1  the  ions  peaked  at  the  appropriate  distance,  but  the 
concentrations  were  far  lower  than  measured.  We  do  not  know  why  the  temperature  effect  did  not  show 
up  here. 


The  difference  between  the  measured  ion  and  neutral  species  used  as  input  data  in  our  modified 
code  runs  compared  with  the  calculated  C3H3+  ion  and  neutral  species  profiles  is  shown  in  Figs.  5  to 
7.  The  calculated  reactive  linear  propargyl  ion,  C3H3+(1)  and  experimental  C3H3+  maximum 
concentrations  are  reasonably  close,  Fig.  5,  but  at  0.5  cm  above  the  burner,  the  calculated  value  is  about 
an  order  of  magnitude  too  high.  It  would  be  more  reasonable  to  compare  the  cyclopropenylium  ion, 
C3H3(c),  with  experiment;  now  the  calculated  value  is  more  than  two  orders  of  magnitude  greater  than 
the  experimental  value  at  e.g.,  0.5  cm.  Most  of  the  neutral  species  are  not  too  different,  Figs.  6  and  7, 
but  two  reactive  neutrals,  C4H2  and  C4H3  are  different  by  about  an  order  of  magnitude  at  0.5  cm,  in 
opposite  directions. 
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In  our  initial  computer  runs,  employing  as  input  data  the  experimental  concentration  profiles 
displayed  in  Figs.  5  to  7,  we  observed  even  less  agreement  with  experiment,  which  we  ascribe  to 
temperature  (equilibrium)  effects,  Fig.  8.  In  the  calculated  ion  profiles,  a  first  peak  appears  far  closer 
to  the  burner  than  is  realistic,  but  second,  smaller,  peaks  are  at  approximately  the  correct  position.  Our 
working  hypothesis  is  that  this  is  a  thermodynamic  temperature  effect.  One  possible  explanation  is  that 
the  ion-molecule  reactions  have  an  activation  energy;  such  reactions  usually  do  not  have  activation 
energies.  If  we  assume  activation  energies  between  80  and  250  kJ/mol  these  primary  peaks  can  be  made 
to  disappear  or  become  fairly  small.  Small  peaks  in  ion  concentration  profiles  have  been  observed  near 
the  burner  experimentally14,15;  the  above  temperature  effect  may  explain  that  dilemma. 

We  ask  why  the  thermodynamics,  which  controls  the  reverse  reaction  rates,  plays  such  an 
important  role.  This  effect  seems  to  be  even  greater  for  neutral  reactions.4,5  We  also  ask  why  the 
calculated  concentrations  are  so  much  smaller  than  the  measured  concentrations.  It  could,  of  course,  be 
that  ion  growth  is  not  the  source  of  larger  ions  in  flames,  but,  we  have  eliminated  many  other 
mechanisms.4,16  Homann  and  associates17  seem  to  be  coming  around  to  recognizing  that  the  source  of 
ions  is  ion  growth  as  we  proposed  years  ago! 18  Maybe  the  experiments  are  off.  It  is  difficult  to  see  how 
they  could  be  off  by  this  amount.  Equivalent  experimental  results  have  been  obtained  by  Homann  and 
associates  in  Germany19,20  and  by  Delfau  and  associates  in  France.21  We  have  not  yet  incorporated  the 
results  of  the  modified  Langevin  theory,  but  we  do  not  anticipate  that  the  changes  will  be  sufficient  to 
bring  the  experimental  and  calculated  concentrations  together. 

Several  rather  "brutal"  modifications  have  been  made  in  an  attempt  to  determine  what  has  to  be 
done  to  obtain  agreement.  We  are  now  examining  these  in  more  detail  to  determine  which  are  reasonable 
and  which  are  nonsense.  Clearly  at  some  point,  if  we  cannot  obtain  agreement  with  reasonable 
assumptions,  we  may  have  to  abandon  the  concept  that  large  ions  grow  from  chemiions.  This  would 
obviously  lead  to  the  conclusion  that  ions  have  nothing  to  do  with  soot  formation.  A  means  of  explaining 
the  concentrations  of  large  ions  will  then  have  to  be  developed.  Because  of  the  present  situation,  we  will 
reexamine  some  of  the  proposed  sources  of  large  ions  which  we  have  previously  reviewed  and  found 
wanting.  From  the  standpoint  of  understanding  the  mechanism  of  soot  formation,  the  removal  of  the 
ionic  mechanism  from  the  arena  does  not  mean  that  current  free  radical  mechanisms  are  to  be  accepted 
by  default.  Closer  examination  of  the  free  radical  mechanisms  seems  to  lead  to  the  same  problems;  the 
models  do  not  agree  with  experiment.4  As  pointed  out  previously  4  the  free  radical  mechanisms  do  not 
predict  the  neutral  species  that  are  observed  in  flames,  and  when  they  do,  they  predict  concentrations 
more  than  an  order  of  magnitude  too  small  for  fairly  small  polycyclic  aromatic  hydrocarbons. 

The  results  of  the  brutal  modifications  to  the  mechanism  were  not  expected,  e.g.,  reducing  the 
ion  recombination  rate  had  no  effect  on  the  profiles!  Further  examination  of  this  set  of  results,  and 
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previous  results,  has  given  us  some  concern  as  to  whether  the  computer  code  is  appropriately  handling 
the  idiosyncratic  requirements  of  ions.  For  example,  the  results  of  the  detailed  computer  modeling  are 
inconsistent  with  conclusions  from  calculating  the  time  to  add  a  carbon  atom  to  the  growing  species4  and 
the  strong  thermodynamic  driving  force  for  the  ionic  mechanism,  see  Section  DC. 

We  are  now  checking  the  code  to  determine  whether  or  not  it  is  appropriately  treating  ions  and 
whether  we  are  using  it  correctly,  and  we  are  working  toward  making  the  above  assumptions  more 
reasonable.  This  will  involve  step  by  step  analysis  of  the  mechanism  to  identify  the  problem  areas.  We 
will  also  evaluate  the  use  of  more  neutral  reactants  in  the  growth  mechanism,  especially  C3H4,  which  in 
some  reaction  steps  has  been  observed  to  be  important.  Toward  simplification,  we  have  based  most  of 
the  mechanism  on  acetylene  because  the  maximum  concentration  of  acetylene  is  about  20  times  greater 
than  the  maximum  concentration  of  diacetylene,  the  next  largest  concentration  of  reactive  species,  and 
700  times  greater  than  the  concentration  of  allene.  When  the  reacting  ion  concentration  is  large,  reaction 
with  C4H2  or  C3H4  in  an  exothermic  reaction  may  be  a  better  path  than  reacting  with  in  an 
endothermic  reaction. 


VIII.  COMPARISON  OF  THE  RATE  OF  CARBON  SPECIES  GROWTH 
FOR  THE  FREE  RADICAL  AND  IONIC  MECHANISMS 


In  the  final  report  on  the  last  contract,1  we  compared  the  rate  of  increase  in  carbon  content  by 
a  free  radical  and  the  ionic  mechanism,  using  experimentally  measured  neutral  and  ion  species,  and  the 
appropriate  rate  coefficients.  That  comparison  has  been  further  tuned  and  has  been  submitted  for 
publication.4  When  the  more  relevant  experimental  data  of  Bittner  and  Howard22  are  used  rather  than 
the  data  of  Vovelle23  (our  burner  and  conditions  were  designed  to  match  those  of  Bittner  and  Howard 
so  such  comparisons  could  be  made)  the  ionic  mechanism  wins;  compare  6.7  /ts  for  the  ionic  mechanism 
with  25  /is  for  the  free  radical  mechanism.  This  is  still  not  persuasive,  recognizing  the  inaccuracies  in 
the  data  and  in  die  rate  coefficients,  but  it  certainly  demands  attention. 

The  experimental  data  on  which  this  analysis  was  made  are  presented  in  Fig.  9.  To  the  extent 
that  Bittner  and  Howard’s  data  can  be  extrapolated  to  beyond  250  u,  the  ion  and  neutral  concentrations 
become  equal  and  the  difference  in  growth  rate  will  be  determined  by:  the  relative  reaction  rate 
coefficients;  the  concentration  of  the  other  reactants-  hydrogen  atoms  and  acetylene  for  the  free  radical 
mechanism,  and  acetylene  for  the  ionic  mechanism;  and  the  number  of  reaction  steps  to  account  for  the 
same  number  of  carbon  atoms  added,  ratio  2/1  for  the  free  radical  mechanism  over  the  ionic  mechanism. 
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It  is  becoming  increasingly  difficult  to  defend 
mechanism. 


a  free  radical  mechanism  against  an  ionic 


IX.  COMPARISON  OF  THE  THERMODYNAMIC  DRIVING  FORCE 
FOR  THE  IONIC  AND  FREE  RADICAL  MECHANISMS 


The  ultimate  driving  force  for  any  sequence  of  chemical  reactions,  such  as  must  occur  for  soot 
formation,  is  thermodynamics.  As  water  flows  downhill,  reactions  proceed  in  the  direction  of  a  negative 
gradient  in  free  energy.  To  compare  this  driving  force  for  ions,  free  radicals,  and  neutral  molecules,  we 
compare  the  free  energy  of  formation  of  these  three  components  of  the  ionic  and  free  radical  mechanism 
of  soot  formation  in  Figs.  10  and  11  for  the  two  temperatures,  1000  and  2000  K. 

Several  features  are  obvious  on  inspection.  First,  the  driving  force  for  ions  is  distinctly  greater 
than  for  free  radicals  or  neutral  molecules.  Second,  the  free  energy  per  carbon  atom  of  either  a  free 
radical,  a  neutral  molecule,  or  an  ion  rapidly  levels  off  as  the  compound  becomes  larger.  Third,  the  ion 
and  neutral  species  free  energies  of  formation  converge  when  the  species  contain  about  30  carbon  atoms, 
i.e.,  a  molecular  weight  of  about  350  to  400  u. 

The  conclusions  to  be  drawn  are  clear.  There  is  very  little  thermodynamic  driving  force  for  the 
growth  of  free  radicals  or  neutral  molecules  compared  to  the  driving  force  for  ions. 


X.  IONS  IN  BENZENE/OXYGEN  FLAMES 


In  a  previous  program,6  we  measured  ion  spectra  of  benzene/oxygen  flames  in  the  same  flames 
and  on  the  same  burner  as  used  by  Bittner  and  Howard22  to  measure  neutral  species  profiles.  This  was 
done  to  have  a  data  base  for  some  day  modeling  that  flame.  We  are  currently  analyzing  that  data  in 
preparation  for  modeling  the  benzene/oxygen  flame.  Some  of  that  analysis  was  presented  at  the  Eastern 
States  Section  Meeting  of  the  Combustion  Institute.24 

Several  classes  of  ions  can  be  distinguished  on  the  basis  of  their  profile  shapes.  C3H3+  is 
assumed  to  be  the  primary  ion.  In  both  nonsooting  ( <f>  =  1.8)  and  sooting  (<j>  =  2.0)  flames,  it  exhibits 
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two  peaks,  a  small  one  near  the  burner  and  a  relatively  large  second  peak  downstream.  It  slowly  decays 
further  from  the  burner.  The  unique  feature  of  the  benzene  flame  is  that  a  large  number  of  oxygenated 
ions,  e.g.,  C6H70+,  are  observed  which  are  not  observed  in  acetylene/oxygen  flames.  Most  of  the  ions 
observed  in  acetylene/oxygen  flames  are  also  observed  in  the  benzene/oxygen  flame. 

The  apparent  stability  of  the  C3H3+  ion  in  the  benzene  flame  contrasts  with  its  behavior  in  the 
C2H2/O2  flame.  In  the  benzene  flame  its  concentration  is  an  order  of  magnitude  greater  than  the 
concentration  of  other  ions;  in  the  sooting  acetylene  flame  it  has  the  same  order  of  magnitude 
concentration.  In  the  benzene  flames  it  decays  more  slowly  than  in  the  acetylene  flame.  The  difference 
can  probably  be  explained  by  the  equilibrium  between  the  linear  reactive  ion  and  the  stable  cyclic 
nonreactive  ion  favoring  the  nonreactive  form  at  the  lower  temperature  in  the  benzene  flame  compared 
to  the  acetylene  flame.  Some  of  the  differences  are  probably  due  to  the  fact  that  the  acetylene  flame  was 
more  fuel  rich,  relative  to  the  equivalence  ratio  for  soot  formation,  than  the  benzene  flames.  Since  the 
total  ion  concentrations  in  the  benzene  flames  are  greater  than  that  in  the  acetylene  flame,  compare  about 
10n  with  7  x  109  ions  cm'3,  the  stability  of  C3H3+  (cyclic)  apparently  does  not  limit  the  formation  of 
large  ions,  but  spreads  their  formation  out  over  a  greater  distance.  These  differences  between  the  two 
flames  will  hopefully  be  explained  by  the  models. 

Recently  Loffler  and  Homann25’26  reported  ion  profile  measurements  in  sooting  benzene/oxygen 
flames.  The  conditions  in  their  experiments  and  in  ours  are  compared  in  Table  I.  In  Table  II  we 
compare  the  ions  observed  in  acetylene/oxygen  and  benzene/oxygen  flames  and  the  ions  observed  at 
AeroChem  and  at  Darmstadt  in  both  flames.  Loffler  and  Homann  concentrate  in  their  paper  on 
oxygenated  ions,  but  report  observed  spectra  in  Loffler’s  thesis.  We  have  thus  used  these  spectra  to 
determine  other  ions  observed  in  their  flames. 

In  any  such  comparison  of  ions,  or  mass  spectrometer  data,  one  should  be  aware  that  the  choice 
of  the  appearance  of  an  ion  contains  a  somewhat  subjective  factor  in  choosing  which  peaks  to  identify 
as  important,  and  an  objective  factor  that  is  dependent  upon  the  sensitivity  and  resolution  of  the  mass 
spectrometer.  The  list  used  here  for  acetylene/oxygen  flames  at  Darmstadt  was  kindly  furnished  by 
Homann28  to  supplement  those  given  in  the  paper.19  Many  other  ions  were  identified  by  Homann;  we 
used  only  those  identified  as  "prominent  peaks" .  Some  of  these  ions,  e.g.,  71, 83,  and  85  u  are  identified 
by  Gerhardt  and  Homann  in  the  acetylene  flame  as  an  oxygen  molecule  attached  to  an  ion.  This  indicates 
reactions  in  the  sampling  cone  of  their  instrument;  these  ions  would  not  be  expected  to  be  stable  in  a 
flame.  Our  experience  has  been  that  the  attachment  of  water,  sometimes  more  than  one  water  molecule, 
to  an  ion  indicates  sampling  cone  problems.  Such  ions  can  be  gotten  rid  of  by  changing  the  cone 
structure,  temperature,  voltage,  or  pressure  behind  the  cone.  Some  of  their  oxygenated  species  may  in 
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fact  be  water  clustered  to  an  ion,  e.g.,  47,  67,  69,  71,  73,  83,  85,  95,  105,  109,  and  119.  This  could 
be  the  case  for  some  of  the  oxygenated  species  observed  in  both  laboratories  in  benzene  flames. 

The  above  consideration  raises  again29  the  sampling  problem  of  ions  which  have  very  high 
reaction  rate  coefficients.  Are  the  identified  ions  the  true  ion  or  an  ion  that  has  been  produced  in  the 
sampling  system,  by  e.g.,  reacting  with  a  species  in  large  concentration,  e.g.,  H20,  as  discussed  above, 
or  H2,  or  H  • ,  or  G,H2?  This  tends  to  muddle  an  already  complicated  problem.  It  means  that  the  ion 
in  the  mechanism  should,  in  some  instances,  be  one  that  can  simply,  in  an  elementary  reaction,  produce 
the  observed  ion.  This  may,  in  fact,  shed  some  light  on  the  problems  we  are  experiencing  with  computer 
modeling  the  mechanism.  Section  VI. 

In  general,  for  the  benzene  flame,  Table  n,  the  agreement  in  observed  ions  between  Loffler  and 
Homann  and  AeroChem  are  consistent,  except  for  C13H9+.  We  observe  this  ion  but  they  do  not.  This 
is  a  very  stable  ion  and  we  both  observe  it  in  large  concentrations  in  acetylene  flames.  Instead  of  this 
ion,  they  observe  C12H90+,  4  u  greater  than  C13H9+,  an  ion  we  observe  only  in  very  small 
concentrations.  Because  we  had  previously  reported  C13H9+  in  benzene  flames,27  Loffler  and  Homann 
were  very  careful  to  confirm  the  fact  that  it  was  not  observed.30  Their  instrumentation  is  better  than 
ours,  but  it  is  difficult  to  see  how  we  could  be  off  in  our  calibration  by  4  u.  If  we  shift  our  mass  scale 
by  4  u  to  force  a  match  at  169  u,  then  our  agreement  with  larger  ions  falls  from  80  ions  to  54  ions. 
Further,  we  would  identify  40  ions  which  they  do  not  observe  and  are  not  observed  in  the  acetylene 
flame.  We  thus  need  to  seek  elsewhere  for  this  difference  in  observation  of  Ci3H9+  at  mass  165. 

One  possibility  is  that  the  two  ions  have  similar  thermodynamic  stabilities,  so  that  the  differences 
in  our  two  experiments.  Table  I,  are  sufficient  to  switch  the  stability  from  one  to  the  other  ion.  To 
determine  the  feasibility  of  this  would  require  estimating  the  thermodynamics  for  169  u,  which  we  will 
do.  This  ion  may  also  be  152  u  with  H20+  replacing  a  H  atom  or  151  u  with  H20+  added,  but  then 
why  does  Darmstadt  not  observe  mass  165? 

The  next  step  in  analyzing  the  benzene/oxygen  flame  will  be  to  calibrate  the  mass  spectrometer 
absolute  sensitivity  by  comparing  the  mass  spectral  data  with  our  Langmuir  probe  data. 
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XI.  PirRIir ATIQNS  AND  PRESENTATIONS 

1.  Keil,  D.G.  and  Calcote,  H.F.,  "Large  Ion  Formation  in  a  Sooting  and  Near  Sooting 
Benzene/Oxygen  Flame,"  Fall  Technical  Meeting,  Eastern  Section:  The  Combustion  Institute, 
Orlando,  FL,  3-5  December  1990. 

2.  Calcote,  H.  F.  and  Gill,  R.  J.,  "Formation  of  Large  Species  in  Sooting  Flames  by  Free  Radical 
and  Ionic  Mechanisms,"  Fall  Technical  Meeting,  Eastern  Section:  The  Combustion  Institute, 
Orlando,  FL,  3-5  December  1990. 

3.  Calcote,  H.F.,  "The  Role  of  Ions  in  Soot  Formation,"  presented  at  International  Workshop, 
Mechanisms  and  Models  of  Soot  Formation,  Ruprecht-Karls-Universitat,  Heidelberg,  Germany, 
29  September-2  October  1991. 

4.  Calcote,  H.F.  and  Gill,  R.G.,  "Comparison  of  the  Ionic  Mechanism  of  Soot  Formation  with  a 
Free  Radical  Mechanism,"  to  be  published  in  Proceedings  International  Workshop,  Mechanisms 
and  Models  of  Soot  Formation,  Ruprecht-Karls-Universitat,  Heidelberg,  Germany,  29  September- 
2  October  1991,  Springer-Verlag. 

5.  Gill,  R.J.  and  Calcote,  H.F.,  "Thermodynamic  Properties  of  Organic  Cations,"  December  1991, 
submitted  to  J.  Phys.  Chem.  Ref.  Data. 


XII.  PROFESSIONAL  PARTICIPATION 

The  Principal  Investigator  on  the  program  is  H.F.  Calcote  who  has  been  responsible  for 
developing  the  reaction  mechanism  and  collecting  the  reaction  kinetics  data.  Dr.  C.H.  Berman  extended 
the  Langevin  theory  to  large  molecular  ions.  Dr.  Fokian  Egolfopoulos  modified  the  Sandia  Flame  Code 
to  handle  ions.  Dr.  R.J.  Gill  has  been  responsible  for  collecting  and  calculating  the  thermodynamics  data 
and  for  running  die  modified  Sandia  Flame  Code.  The  benzene  flame  experimental  data  were  previously 
collected  by  Dr.  D.G.  Keil  who  has  assisted  in  reducing  them  in  preparation  for  modeling.  Leslie  Van 
Hoose  and  I.E.  Palmer  have  assisted  with  detailed  calculations  and  Helen  Rothschild  has  assisted  with 
graphics  and  editing.  We  all  acknowledge  fruitful  discussions  with  Dr.  W.  Felder. 
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XIV.  NOMENCLATURE 


E 

E„ 

h 

k 

n 

“v 

% 

r 

rc 


rmin 


r 


* 


x 


y 


radius  of  a  disk 

hard  sphere  neutral  radius 

radius  of  a  sphere 

impact  parameter,  minimum  center  to  center  distance  between  two  molecules  if  field 
forces  did  not  affect  their  trajectories 
electric  field  strength 

electric  field  strength  at  the  location  of  the  neutral  reactant 
thickness  of  a  disk 
collision  rate  coefficient 
total  number  density 

number  of  neutral  reactant  mole  molecules  having  a  given  velocity  v 
electric  charge  on  ion 
point  electrical  charge 

distance  measured  from  the  center  of  the  ion 

critical  radius,  r,  in  Langevin  theory  for  which  neutral  species  describes  an  orbital  motion 
around  the  ion.  Collision  occurs  when  r  <  rc 

minimum  separation  distance  between  ion  and  neutral,  measured  from  the  center  of  the 
ion 

separation  distance  from  center  of  ion  to  center  of  molecule  at  which  the  electric  field 
strength  approaches  infinity 
mean  molecular  velocity 

molecular  velocity  of  neutral  reactant  uninfluenced  by  the  ion 
distance  from  center  of  sphere  to  point  where  field  is  measured 
distance  from  center  of  sphere  to  the  external  point  charge 


7  =  polarizability 

7j  =  polarizability  of  ion 

7„  =  polarizability  of  neutral 

0  =  angular  rate  for  neutral  molecule  orbiting  ion 

/x  =  reduced  mass 

0  =  angular  rate  for  neutral  molecule  orbiting  ion 
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TABLE  I 


COMPARISON  OF  BENZENE/OXYGEN  FLAMES  OBSERVED 
AT  AEROCHEM  AND  AT  DARMSTADT 

AeroChem  Darmstadt 


Equivalence  Ratio  1.8,  2.0 

+  30  mol  96  argon 

Unburned  Gas  Velocity,  cm  s-1  50 

Pressure,  kPa  2.67 

Threshold  Soot  Formation,  <t>  1.90 


1.9,  2.0 

42 

2.67 

1.83 


Temperature  Peak,  K: 
at  mm  above  burner: 
for  equivalence  ratio: 


1900 

12.9 

1.80 


2094 

11.6 

1.95 
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TABLE  n 

COMPARISON  OF  IONS  OBSERVED  IN  ACETYLENE/OXYGEN  AND 
BENZENE/OXYGEN  FLAMES  AND  IONS  OBSERVED  AT  AEROCHEM  AND  AT 
HOMANN’S  LABORATORY  IN  DARMSTADT,  GERMANY 


ION 

ACETYLENE/OXYGEN 

BENZENE/OXYGEN 

MASS 

ION  FORMULA 

AERO* 

DARMS0 

AERO® 

DARMS® 

19 

H30+ 

X 

24 

Q>  + 

X 

X 

31 

ch3o+ 

X 

33 

ch5o+ 

X 

39 

c3h3+ 

X 

X 

X 

X 

41 

C2HO 

X 

X 

43 

C2H30+ 

X 

X 

X 

47 

C2H7o+ 

X 

49 

c4h+ 

X 

51 

C4H3  + 

X 

X 

X 

53 

c4h5+ 

X 

X 

X 

X 

55 

c3h3o+ 

X 

57 

C4H9\  C2HtO+ 

X 

X 

63 

c5h3+ 

X 

X 

X 

X 

65 

c5h5+ 

X 

X 

67 

c4k3o+ 

X 

X 

X 

69 

c4h5o+ 

X 

71 

c3h3+-  Oj 

X 

73 

c6h+,  c4h9o+ 

X 

X 

75 

c6h3+ 

X 

X 

X 

77 

c6h5+ 

X 

X 

X 

79 

c6h7+ 

X 

X 

X 

81 

c5h5o+ 

X 

X 

X 

83 

c4h3+-  O2 

X 

85 

C4H5+-  O2 

X 

87 

c7h3+ 

X 

X 

X 

89 

c7h5+ 

X 

X 

X 

91 

c7h7+ 

X 

X 

X 

X 

95 

c6h7o4 

X 

X 

X 
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TABLE  II  (continued) 


ION 

ACETYLENE/OXYGEN 

benzene/oxygen 

MASS 

ION  FORMULA 

AERO4 

DA  RMS0 

AERO0 

DARMSa 

97 

C8H+ 

X 

X 

99 

CgH3  + 

X 

X 

X 

103 

c8h7+ 

X 

X 

X 

X 

105 

c7h5o+ 

X 

X 

X 

109 

c7h9o+ 

X 

X 

X 

111 

c9h3+ 

X 

X 

X 

115 

c9h7+ 

X 

X 

X 

X 

117 

c9h9+ 

X 

X 

119 

c8h7o+ 

X 

X 

X 

121 

c8h9o+ 

X 

X 

123 

C8HuO+ 

X 

X 

X 

124 

c8h12o+ 

X 

X 

126 

C,oH6+ 

X 

X 

129 

c10h9+ 

X 

X 

X 

X 

131 

C9H70+,  C10Hn  + 

X 

X 

X 

133 

c<>h9o+ 

X 

X 

135 

cuh3+ 

X 

139 

c„h7+ 

X 

X 

141 

c„h9+ 

X 

X 

X 

X 

143 

c10h7o+ 

X 

X 

X 

145 

c10h9o+ 

X 

X 

X 

147 

c12h3+ 

X 

152 

c12h8+ 

X 

X 

153 

c12h9+ 

X 

X 

X 

155 

c„h7o+,  c12hu+ 

X 

X 

157 

cuh9o+ 

X 

X 

159 

c13h3+ 

X 

165 

c,3h9+ 

X 

X 

X 

169 

c  I2^9o0  + 

X 

170 

Cl2H10O+ 

X 

X 

171 

Ci4H3+ 

X 

177 

c14h9+ 

X 

X 

X 

179 

C,4Hu  + 

X 

X 

181 

c13h9o+ 

X 

X 

183 

c15h3+ 

X 
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TABLE  II  (continued) 


ION 

ACETYLENE/OXYGEN 

BENZENE/OXYGEN 

MASS 

ION  FORMULA 

AERO® 

DARMS0 

AERO* 

DARMSa 

189 

c15h9+ 

X 

X 

X 

191 

Ci5Hn  + 

X 

X 

X 

192 

c14h8o+ 

X 

X 

193 

c14h9o+ 

X 

X 

203 

C,6Hn  + 

X 

X 

X 

204 

Cl5H80+ 

X 

X 

205 

c15h9o+ 

X 

X 

215 

Ci7Hn  + 

X 

X 

X 

X 

217 

c16h9o+ 

X 

X 

219 

c,6huo+ 

X 

X 

227 

C,8Hn  + 

X 

X 

X 

X 

228 

c17h8o+ 

X 

X 

231 

n 

X 

o 

+ 

X 

X 

239 

Ci9HU  + 

X 

X 

X 

X 

240 

c18h8o+ 

X 

X 

241 

c18h9o+ 

X 

243 

n 

CO 

X 

o 

+ 

X 

X 

250 

C2oHio+ 

X 

251 

C2oH„  + 

X 

X 

X 

253 

c19h9o+ 

X 

X 

255 

c19huo+ 

X 

263 

C2lH„  + 

X 

X 

X 

X 

265 

C2oH90+,  02^,3+ 

X 

X 

267 

C2oHnO+,  C21Hl5+ 

X 

276 

C22H12  + 

X 

277 

C22H13  + 

X 

X 

X 

279 

c21huo+,  c22h15+ 

X 

X 

287 

C23Hu  + 

X 

289 

C23H13  + 

X 

X 

X 

293 

C23H,30+ 

X 

X 

300 

C24H12+ 

X 

301 

C24H13  + 

X 

X 

X 

X 

303 

C23H„0+ 

X 

X 

313 

C25H13  + 

X 

X 

X 

317 

C24h,30+ 

X 
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TABLE  II  (continued) 


ION 

MASS 

ION  FORMULA 

ACETYLENF/OXYGEN 
AERO*  DARMSb 

BENZENE/OXYGEN 
AEROc  HARMS0 

324 

C26HI2  + 

X 

325 

C26H13  + 

X 

X 

X 

329 

C25H130+ 

X 

337 

C27H13+ 

X 

X 

X 

341 

C,6Hi30+ 

X 

348 

C28H,2+ 

X 

349 

C28H13  + 

X 

350 

C28Hi4+ 

X 

351 

C28H15  + 

X 

X 

X 

353 

C27h13o+ 

X 

361 

C29H,3+ 

X 

X 

X 

X 

363 

C28HuO+ 

X 

X 

X 

365 

c,8h13o+ 

X 

X 

374 

c30h14+ 

X 

375 

c30h13+ 

X 

X 

377 

C29H130+ 

X 

X 

385 

c3,h13+ 

X 

386 

C3iH14+ 

X 

X 

387 

c3,h15+ 

X 

X 

391 

c30h15o+ 

X 

398 

c32h14+ 

X 

399 

c32h15+ 

X 

X 

X 

403 

c31h15o+ 

X 

411 

c33h15+ 

X 

X 

X 

X 

415 

c32h13o+ 

X 

422 

cwh14+ 

X 

423 

c34h15+ 

X 

X 

X 

X 

427 

c33h15o+ 

X 

435 

c33h13+ 

X 

X 

X 

X 

439 

C34Hi30+ 

X 

446 

C^H  14+ 

X 

447 

C36H15+ 

X 

X 

X 

X 

459 

c37h15+ 

X 

X 

X 

X 

470 

c38h14+ 

X 

471 

c38h13+ 

X 

X 

X 
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TABLE  II  (continued) 


ION 

ACETYLENE/OXYOF.N 

BENZENE/OXYGEN 

MASS 

ION  FORMULA 

AERO4 

DARMS0 

AERO* 

DARMS" 

472 

C38H,6+ 

X 

473 

C38H17  + 

X 

X 

483 

C39^15  + 

X 

X 

X 

X 

497 

QoHl7+ 

X 

502 

C4iH15+ 

X 

506 

C4,H15+ 

X 

X 

507 

C4iH15+ 

X 

508 

c41h16+ 

X 

X 

509 

c41h17+ 

X 

X 

X 

520 

c42h16+ 

X 

521 

c42h17+ 

X 

X 

X 

533 

C43H17  + 

X 

X 

X 

X 

544 

X 

545 

C44H,7+ 

X 

X 

X 

X 

557 

C45Hi7+ 

X 

X 

X 

X 

568 

C4«Hi6+ 

X 

569 

Q6Hl7+ 

X 

581 

C47H17+ 

X 

“AeroChem,  Ref.  27. 
bDarmstadt,  Refs.  19,  28. 
°This  work. 
dRefs.  25,  26. 
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El-  hard  sphere 


langevtn  collision rate^ 

ASS™«  ******  “  ”«««l  reac4, 
and  a  single  charge  on  sDhere 


ssumes  acetylene  as  neutral  reactant 
and  a  single  charge  on  sphere 
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ION  DIAMETER,  nm 


FIGURE  3  EFFECT  OF  ION  DIAMETER  ON  THE  RATIO  OF 
THE  EXPANDED  LANGEVIN  TO  THE  HARD  SPHERE  COLLISION  RATE 
Assumes  acetylene  as  the  neutral  reactant 
and  a  single  charge  on  sphere 


DISTANCE,  cm 


FIGURE  4  COMPARISON  OF  EXPERIMENTAL  WITH  CALCULATED  ION 
PROFILES  FROM  LAST  RUN  BY  FRENKLACH  AND  WANG 
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DISTANCE,  cm 

FIGURE  5  COMPARISON  OF  EXPERIMENTAL  SPECIES  PROFILES  USED  AS 
INPUT  DATA  TO  THE  MODIFIED  FLAME  CODE  WITH  CALCULATED 
VALUES  FROM  LAST  RUN  BY  FRENKLACH  AND  WANG 
- Calculated; - Experimental 


DISTANCE,  cm 


FIGURE  6  COMPARISON  OF  EXPERIMENTAL  SPECIES  PROFILES  USED 
AS  INPUT  DATA  TO  THE  MODIFIED  FLAME  CODE  WITH  CALCULATED 
VALUES  FROM  LAST  RUN  BY  FRENKLACH  AND  WANG 
- Calculated; - Experimental 
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DISTANCE,  cm 

1GURE  7  COMPARISON  OF  EXPERIMENTAL  SPECIES  PROFILES  USED 
S  INPUT  DATA  TO  THE  MODIFIED  FLAME  CODE  WITH  CALCULATED 
VALUES  FROM  LAST  RUN  BY  FRENKLACH  AND  WANG 
- Calculated; - Experimental 


DISTANCE,  cm 

FIGURE  8  TYPICAL  CALCULATED  ION  PROFILES  USING  THE 
MODIFIED  FLAME  CODE 
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ENERGY  OF  FORMATION 

kJ/mol/atom  Q  MOLE  FRACTION 


NUMBER  CARBON  ATOMS 

FIGURE  10  FREE  ENERGY  OF  FORMATION  PER  CARBON  ATOM  FOR  IONS, 
FREE  RADICALS  AND  MOLECULES  AT  1000  K 
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FIGURE  1 1  FREE  ENERGY  OF  FORMATION  PER  CARBON  ATOM  FOR  IONS, 
FREE  RADICALS  AND  MOLECULES  AT  2000  K 
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